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Abstract 

We consider the analogue of the 4th Appelrot class of motions of the Kowalevski top for 
the case of two constant force fields. The trajectories of this family fill the four-dimensional 
surface D in the six-dimensional phase space. The constants of three first integrals in in- 
volution restricted to this surface fill one of the sheets of the bifurcation diagram in R 3 . 
We point out the pair of partial integrals to obtain the explicit parametric equations of this 
sheet. The induced system on D is shown to be Hamiltonian with two degrees of freedom 
■ having the thin set of points where the induced symplectic structure degenerates. The region 

of existence of motions in terms of the integral constants is found. We provide the separation 
of variables on O and the algebraic formulae for the initial phase variables. 
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1 Preliminaries 

The equations of motion of the Kowalevski top in two constant fields, expressed in the reference 
system Oeie 2 e 3 of the principal axes of inertia at the fixed point O, 



2u)i = UJ 2 UJ 3 + /?3, 2^2 = -OJ1CO3 - «3, V3 = OC2 — Pl, 

di = a 2 u 3 - a 3 u 2 , $1 = (3 2 uj 3 - (3 3 uj 2 , 

a 2 = a 3 u x - aicj 3 , $ 2 = #3^1 - 

a 3 = o>iuj 2 — a 2 uJi, $ 3 = (3\uj 2 — (3 2 uji 



'l.D 



can without loss of generality be restricted to the phase space P 6 C R 9 (w, a, (3) denned by the 
geometric integrals 

a\ + a\ + a\ = a 2 , 0* + f3 2 2 + /3f = b 2 , a t Px + a 2 fc + a 3 fa = (1.2) 
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(see [T] for details). This system is completely integrable due to the existence of three integrals in 
involution [21 [3] 

H = u)\ + uj\ + -uj - (a l + (3 2 ), 



K = (uj - u 2 2 + oti - /3 2 ) 2 + (2c^ 2 + a 2 + A) 2 , (1.3) 



G = -{Ml + M 2 ) + -cu 3 M 7 - - a 2 p 2 , 



where 



M a = 2u)\a,\ + 2u; 2 a 2 + ^3^3, 

Mp = + 2u 2 p 2 + u 3 p 3 , (1.4) 

M 7 = 2ui(a 2 fc - a 3 p 2 ) + 2^ 2 (a 3 /3i - a.xPz) + u 3 {a x p 2 - a 2 Pi)- 

For the general case 

a>b>0 (1.5) 

the explicit integration has not been found yet. It is natural to study first the invariant subman- 
ifolds in P 6 such that the induced system has only two degrees of freedom. It is proved in [U Q] 
that there exist only three submanifolds 9Jt, 91, of this type. The union Wl U 01 U D coincides 
with the set of critical points of the integral map 

H x K x G : P 6 -> R 3 . (1.6) 

If 6 = (the classical Kowalevski case [5]), then the critical set of the map (11.61) consists of 
the motions that belong to the so-called four Appelrot classes [6]. The set 9JI, first found in [2] 
as the zero level of the integral K, generalizes the 1st Appelrot class. The phase topology of the 
system induced on DJl was studied in [7]. The dynamical system on 01 generalizing the 2nd and 
3rd Appelrot classes was explicitly integrated in [8]. The present work considers the restriction of 
the system (II. ip to the invariant subset O. 

Introduce the complex phase variables [9] (i 2 = — 1): 



xi = («i - p 2 ) + i(a 2 + /3i), x 2 = (ai - f3 2 ) - i(a 2 + Pi), 
yi = (aa + /3 2 ) + «(a 2 - Pi), y 2 = (a x + (3 2 ) - i(a 2 - Pi), 
z 1 =a 3 + ip 3 , z 2 = a 3 - z/3 3 , 
iwj = co>i + za; 2 , u> 2 = uj 1 — iuj 2 , w 3 = tu 3 . 



;i-7) 



The system (II .ip takes the form 

x[ = —XlWz + ZiWi, x' 2 = x 2 w 3 - Z 2 W 2 , 

V'l = -ViW 3 + z 2 w u y' 2 = y 2 w 3 - z x w 2 , , lg , 
2z[ = x 1 w 2 - y 2 w 1 , 2z' 2 = -x 2 w 1 + yiw 2 , 
2w[ = -(wiw 3 + zi), 2w' 2 =w 2 w 3 + z 2 , 2w' 3 = y 2 -y 1 . 

Here the prime stands for d/d(it). 

The set 0, by the definition given in the work pQ, includes as a proper subset the following 
points 

wi = w 2 = 0, zi = z 2 = 0. (1.9) 
The invariant relations (I1.9P lead to the family of pendulum motions first found in [I] 

a = a(ex cos 9 — e 2 sin 0) , /3 = ± fe(ei sin + e 2 cos #) , 

p, i l d» ^ , , - • fl (L10) 

a x p = ± aoe 3 , u? = — ea, -— = — (a ± o )sm9. 

dt ' dt 2 K 1 
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The corresponding values of the integrals (I1.3P satisfy one of the following 

g = ±abh, = (a ^ 6) 2 , h^-(a±b). (1.11) 

In the sequel we use functions and expressions having singularities at the points (11.91) . There- 
fore, by default, we exclude the trajectories (11.101) . The remaining part of O can be described by 
the pair of equations 

#i = 0, R 2 = 0, (1.12) 

where 

w 2 x 1 + wiy 2 + w 3 z 1 wix 2 + w 2 yi + w 3 z 2 



Ri 



Wi w 2 
w 2 z\ W\Z\ 



R 2 = (w 2 z 1 + w 1 z 2 )w 3 + 1 h wiw 2 {yi + y 2 > + xiw 2 + x 2 w 1 w 3 + (1-13) 

L wi w 2 

w\x\Z\ w\x 2 z 2 . 

H 1 h X 1 Z 2 W 2 + X 2 2;iWi + (W7iZ 2 - W2Zi)[yi ~ Vl)- 

Wi w 2 

For the derivatives we have, in virtue of (11.81) . 

R[ = k 2 R 2 , R' 2 = k>iRi, 

Ki = [(WiW 2 W 3 + Z 2 W 1 + Z X W 2 ) 2 + W 1 W 2 (x 2 w\ + XiwD] , K 2 



;i.l4) 



2w 1 w 2 ~ 2w\W 2 



The equations (11. 14j) straightforwardly prove that the set (jl,12p is preserved by the phase 
flow (11.11) . Besides, it follows that the Poisson bracket {Ri,R 2 } is a partial integral on O. The 
subset in O defined by the equation 

{R u R 2 } = (1.15) 

is the set of points at which the induced symplectic structure is degenerate. Below we calculate 
{Ri,R 2 } using the appropriate integrals on O and show that the set (11.151) is of measure zero. 
Thus, the dynamical system induced on D by the system (11.11) is almost everywhere a Hamiltonian 
system with two degrees of freedom. In particular, almost all its integral manifolds consist of two- 
dimensional Liouville tori. 



2 Partial integrals 

Recall that in the classical Kowalevski case {(3 = 0) there exists the momentum integral 

L=hu;.cx (I = diag{2,2,l}). (2.1) 

Then the integral G becomes equal to L 2 . Let i denote the constant of the integral (12.10 . In 
Appelrot's notation the 4th class of especially remarkable motions is defined by the following 
conditions. 

(i) The second polynomial of Kowalevski has a multiple root. One of the Kowalevski variables 
remains constant and equal to the multiple root s of the corresponding Euler resolvent (p(s) = 
s(s-h) 2 + (a 2 -k)s-2£ 2 : 

<p(s) = 0, cp'{s) = 0. (2.2) 
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(ii) Two equatorial components of the angular velocity are constant: u% = —i/s, uj 2 = 0. 
Given f3 = 0, this fact can be written in the form 

\lo • Ct 

= -s, Iu>-/3 = 0, Iu>-e 2 = 0. (2.3 

Io;-ei 

The next statement establishes the conditions similar to (I2.3P for the generalized top. 
Theorem 1. On each trajectory belonging to D the ratios 

Ioj-ei' loj-e 2 

are constant and equal to each other. 
Proof. Let 

M = Iu>, Mj = luf-ej (j = 1, 2, 3). (2.4) 
Due to (11. 4p we also have M a = luj-ct, Mg = 1uj-(3. Then the first equation ( 11. 121) yields 

M a Mg . . 

—rj- — —j- = 0. 2.5 

Mi M 2 K J 

Introduce the function 

M n .\l. ■ .U..U, 
Ml + Mf 

and calculate its time derivative: 

dS {Ml + M 2 )u 3 + 4a 3 Mi + 4/3 3 M 2 



dt ' 2{Ml + Ml) 2 



(M a M 2 — MgMi). 



In virtue of (12.51) the right hand side is identically zero. Therefore S is a partial integral on O. 
Denote by s the corresponding constant: 

M a M 1 + M g M 2 

WTmI = ~ s - (2 - 6) 

From (jZS), (EH) we obtain 

M a = -sM u Mg = -sM 2 (2.7) 

with the constant value s. □ 

Note that according to (12 .7p the function 5* can be written in either of the representations 

s= 1 / Mg + jMg ^ M a — iMg \ _ ! M a + lMg = ! M a -tMg 

4 \ ui + iuj 2 iO\ — / 2 cg>i + zcj 2 2 io\ — iu 2 
Theorem 2. On the set D the system (II. ip /ias t/ie partial integral 

T = ^(M a Mi + MgM 2 ) - 2(a l [3 2 - a 2 ft) + a 2 + b 2 . (2.9) 



Proof. The time derivative of T 



dT 1 , 

— = -u 3 (M a M 2 - M P M X ) 



vanishes on O due to (12.51) . □ 
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Denote by r the constant of the integral T. 

Remark 1. In the work [T| equations similar to fll.l2f) were derived from the condition that the 
function 

2G+(r -a 2 -b 2 )H + sK (2.10) 

with Lagrange's multipliers s,t has a critical point on P 6 . Using the coordinates (11.71) we have 
from (23), (ED 

_ Wi(x 2 Wi + yiw 2 + ^ 2 w 3 ) + w 2 (y 2 w 1 + x x w 2 + Z1W3) 
x 2 wi + yiw 2 + -22^3 y 2 wi + x\w 2 + ziw 3 



2w 2 2wi 

T = hw^X^! + yxW 2 + Z 2 W 3 ) + W 2 (y 2 W! + XiW 2 + Z!W 3 )} + X!X 2 + Z!Z 2 = 

2 (2.12) 
= -2S w x w 2 + xix 2 + z x z 2 . 

Comparing (12. lip . (12 . 1 2[) with the expressions for s,t given in pQ we see that on O Lagrange's 
multipliers in (12.101) coincide with the constants of the integrals S, T introduced here. 

Supposing ( 11 .51) . introduce parameters p, r (p > r > 0) such that 

p i = a 2 + b 2 , r 2 = a 2 -b 2 . (2.13) 

Let h, k,g denote the constants of the general integrals ( 11.31) . Then according to Remark (TJ from 
the results of the work [Ij we obtain that on D 

p 2 — r , r 2 — 2p 2 r + r 4 p 4 — r 4 1 , 9 

h= — +a ' k= — h — +T - 9 = — + 2 (p - T)a - (2 - 14) 

These equations can be considered as parametric equations of the corresponding bifurcation sheet 
of the integral map (II .61) . Eliminating r we have 

#0 = 0, V'(s) = 0, (2.15) 

where 

4 4 

V»(s) = s 2 {s - h) 2 + (p 2 - k)s 2 - 2gs + P-ZJL. 

If /3 = (p 2 = r 2 = a 2 ), then if)(s) = sip(s) and the conditions (12. 15[) turn to (I2.2p . Therefore, 
the set of trajectories belonging to D is the generalization of the set of the especially remarkable 
motions of the 4th Appelrot class. 

3 Parametric equations of integral manifolds 

Due to the equations (12.141) the functions S, T form a complete system of first integrals on D. In 
particular, the equations of the integral manifold 

{CeP 6 : H(Q = h,K(Q = k,G(Q = g} 

in this class of motions are equivalent to the relations (11.121) and the equations 

S = s, T = t. (3.1) 
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Using (12.51) . (12. lip we replace the equations R\ = 0, S = s by 

(y 2 + 2s)w 1 + x x w 2 + z x w 3 = 0, , 2 . 
x 2 w x + (yi + 2s)w 2 + z 2 w 3 = 0. 

At the same time, from fl 1 . 1 3 [) and (12.121) the system R 2 = 0, T = r is equivalent to 

x 2 ziwi + xiz 2 w 2 + (r - xix 2 )w 3 = 0, (3.3) 

2s wiw 2 - (x x x 2 + ziz 2 ) + r = 0. (3.4) 

To obtain a closed system of equations we must add the geometric integrals (jl.2p . In terms of the 
variables (11. 7p we get 

+ xij/2 = ^2 + x 2Vi = r 2 , (3.5) 

^1^2 + ym + 2z 1 z 2 = 2p 2 . (3.6) 

The variables (11.71) by definition satisfy 

x 2 = xi, y 2 = yi, z 2 = zi, w 2 = wi, w 3 e R, (3.7) 

thus forming a space of real dimension 9. Seven relations (13. 2 j) - (13.61) define the integral manifold. 
In the case of independency of the integrals S, T this manifold is two-dimensional and consists of 
Liouville tori filled with quasi-periodic motions. 
Let 

x = \/xix 2 , z = a/ ziz 2 , i = 2swiw 2 . (3.8) 

Then from (E3D 

£ = x 2 + 2 2 - r. (3.9) 

(r 2 - Xiy 2 )(^ 2 - Z22/1) = (3-10) 

yi y 2 = 2p 2 -x 2 -2z 2 . (3.11) 

r 2 (xi?/2 + S23/1) = r 4 + 2pV - (x 2 + z 2 ) 2 . (3.12) 



From (I3.5p . (13. 8p we obtain 
and the equation (I3.6P yields 
Hence 

Rewrite (I3.5P in the form 



[zi + z 2 ) 2 = 2r 2 - (xiy 2 + x 2 y x ) + 2z 2 , 
[zi - z 2 f = 2r 2 - (xi?/2 + x 2 yi) - 2z 2 



(3.13) 



and substitute (13.121) to obtain 

r 2 { Zl + z 2 ) 2 = $!, r\z x - z 2 f = $ 2 , (3.14) 

where 

$ x = ( x 2 + z 2 + r 2 ) 2 - 2(p 2 + r 2 )x 2 = (f + r + r 2 ) 2 - 2(p 2 + r 2 )x 2 , . , 

$ 2 = ( x 2 + z 2 - r 2 ) 2 - 2{p 2 - r 2 )x 2 = (£ + r - r 2 ) 2 - 2(p 2 - r 2 )x 2 . ^ ^ 

Note that the equilibria = of the system (11.11) are included in the family of motions 
(11.101) . On the rest of the trajectories in D the determinant of the equations (13.21) . (I3.3P in Wj 
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(j = 1,2,3) vanishes identically. Calculate this determinant and eliminate zf,z 2 from (13.51) and 
the product y±y 2 from (13.111) to obtain 



2s[(r 2 x 1 - ryi) + (r 2 x 2 - ry 2 )] = -r 2 {x x y 2 + x 2 yi) + 
+2[2s 2 (r - x 2 ) + p 2 (r + x 2 ) - r(x 2 + z 2 )}. 



(3.16) 



On the other hand, the direct calculation in view of (I3.8p . (13.111) gives 

(r 2 x 1 - ry 1 )(r 2 x 2 - ry 2 ) = r A x 2 + r(2p 2 - x 2 - 2z 2 ) - r 2 r(x l y 2 + x 2 y\). (3.17) 

Denote 

a = r 2 - 2p 2 r + r 4 , x = ^ > 0. (3.18) 
From the second relation (12.141) we have the identity 

4s\ 2 = a + 4s 2 r. (3.19) 

Introduce the complex conjugate variables 

III = r 2 x x - ryi, \i 2 = r 2 x 2 - ry 2 . (3.20) 

Eliminating X\y 2 + x 2 y\ from (13. 16j) . (13.171) with the help of ( 13.121) we come to the system 



2s(fi 1 +fi 2 ) = e-^ 2 (x 2 -r)-a, 
\i\\x 2 = r£ 2 + ax 2 — ra. 

Choose 



(3.21) 



n\ = y/2snx - 4s 2 r, /j,* = ^2sfx 2 - As 2 r (3.22) 
to be complex conjugate. Then the system ( 13.211) takes the form 

+ A*2) 2 = ^i, (/i!-/4) 2 = * 2 , (3.23) 

where 

= £ 2 - As 2 {x + x) 2 , ^2 = ? -As 2 (x-x) 2 - (3.24) 

As was to be expected from dimensional reasoning, all phase variables could be expressed in 
terms of two almost everywhere independent auxiliary variables. The above formulae emphasize 
the special role of the pair (x, £). 

From (13T22|) . (13T23D we find 

/ii = 2sr + ^( v / ^+ v 7 ^) 2 , fi 2 = 2sr + ^{^/W 1 - ^T 2 ) 2 . (3.25) 
8s 8s 

Further on the calculation sequence is as follows. From (13.141) we find zi, z 2 . Multiplying the 
equations (I3.20p by x 2 , x\ respectively and using (13.51) we get 

x 2 fii = r 2 x 2 — rx 2 yi = r 2 {x 2 — r) + tz 2 , 
Xifi 2 = r 2 x 2 — TX\y 2 = r 2 {x 2 — r) + rzf, 
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whence X\, x 2 are found. Substituting these values back to (I3.20p we find yi, y 2 . As a result, after 
some obvious transformations we obtain the following expressions for the configuration variables: 



X 2 



2s 4r 4 (x 2 - r) + r( + 



r 2 16 S 2 T + ( V /^ 1 - 
2s4r 4 (x 2 - r) + r(y^ 



^^2 



(3.26) 



2s 



16s 2 r + ( v /* 1 
4[2r£ - t(x 2 - r) - 



^T)2 



2/2 = 2s 



16s 2 r + 
4[2t£ - r(x 2 - r) + a] 



*2_) 2 
$1 ■+ 



2l 



2r vv 1 



16s 2 r + 



" 2 = 2^ ( 



^ 2 ) 2 



(3.27) 



(3.28) 



Note that all radicals are algebraic. The formal choice of signs in ( 13.26l) -( |3.28l) is defined by 
the initial choice in the expressions for fj,i, /i 2 , z±, z 2 . The polynomials tyj (j = 1,2) obviously 
split to multipliers linear with respect to x, £. Then, typically, the projection of an integral 
manifold onto the (x,£)-plane has the form of a quadrangle. Fix any inner point (x, £) of such 
projection. Then the expressions (I3.26I) - (I3.28I) define eight points of the configuration space with 
different set of signs of the radicals v$2, a/^V^2- 

To find W3, use the energy integral. On account of (I2.14p it takes the form 

2s w 2 + 4s wiw 2 - 2s(y 1 + y 2 ) = 4s 2 + 2p 2 - 2r. 
Substitute 2swiW 2 from (13.41) and replace 2p 2 by its expression from (13.61) to obtain 

2s w 2 = D, 

where 



(3.29) 

D = ( Vl + 2s) (y 2 + 2s) - x x x 2 (3.30) 

is the determinant of (13.21) with respect to wi,w 2 . In particular, this determinant vanishes along 
with W3. Due to this fact we can express wi,w 2 either as linear functions of w 3 , or in inverse 
proportion to w 3 : 



w 2 



x\z 2 - (yt + 2s) Zi 

(yi + 2s)(y 2 + 2s)-x 2 
x 2 z x - {y 2 + 2s) z 2 

( yi + 2s)(y 2 + 2s)-xi 



w 3 



w 3 



x\z 2 - (yi + 2s)z 1 

2s w 3 
x 2 z x - (y 2 + 2s)z 2 

2s W3 



(3.31) 



With two possibilities of the sign choice for 1^3 in (13. 29f) we have that the inverse image in 
the integral manifold {5* = s, T = r} H O of a generic point (x, £) contains 16 points. 
We also need the explicit formula for w 3 in terms of x, £. From (13.271) we write 



y x + 2s = 2s 
y 2 + 2s = 2s 



4[2r^-r(x 2 -r)+4s 2 x 21 



16s 2 r + (V^! - 
4[2r£ - t(x 2 - t) + 4s 2 x 2 ] - (V^i + 



^2? 
^2? + 



#1 + 



^ 2 



(3.32) 



16s 2 r + (y/^! + y/V 2 ) 
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Note that 

[16s 2 r + ( - \/W 2 ) 2 } [16s 2 r + (yf$l + \/W 2 ) 2 } = 64s 2 (r£ 2 + ax 2 - ra). (3.33) 

Then 



= V /$ 1 $ 2 ^ 1 ^ 2 -P 
2(r^ 2 + (xx 2 - ra)' 



where 



P( X; £) = ^ + 2t£ 3 + 2[(t-2s 2 -p 2 )x 2 -t(2s 2 -p 2 )-r 4 }£ 2 - 

-8s 2 [(r - 2 X 2 )x 2 + T X 2 }i - 4s 2 (x 2 - X 2 )[2(r - p 2 )x 2 - (r 2 - r 4 )]. 

Let 

Q(x, = ( f + r + 2s 2 - p 2 ) 2 - 4aV + r 4 - (2s 2 - p 2 ) 2 , 

P 2 (x, = - 2xQ(x, £) v 7 ^ 2 + ro 2 " ra. 

Direct calculation proves the identity P1P2 = $i$2^ / i x I / 2- On the other hand, Pi + P 2 = 2P. 
Therefore, 

D (v^-v 7 ^) 2 



4(r£ 2 + ax 2 - ra) ' 

whence 

w 3 = i — ; ^ (3.34) 

2^25 {re + ax 2 - ra) 

Below we use one more representation of u>i, u>2 not depending on io 3 . The system (13.21) yields 

^2^1 = - (?/ 2 + 2s)wiw 2 , x\w\ = -z 1 w 2 w 3 - (yi + 2s)w 1 w 2 . 

Substitute the expressions for WiW 3 , w 2 w 3 found from the inverse proportion in (I3.3ip and eliminate 
WiW 2 by (13 .4p to obtain 

2sx 2 w\ = — (jUi — 2st) — 2sx 2 , 2sxiwl = —(fj> 2 — 2sr) — 2sa; 2 . 

Then from (I3.25P we find 

w 1 = -^—{^Q[+^Q 2 ) 1 w 2 = -^—{^Q' 1 -^Q 2 ), (3.35) 

where 

©i(z, = (£ - 2sa;) 2 - 4s 2 X 2 , e 2 (x, = (£ + 2sa;) 2 - 4s 2 X 2 . (3.36) 

Note that the products ©162 and ^1^2 coincide. The choice of the signs in (13.35j) is deter- 
mined by the condition (y/x^Wi)(^Jx~iW 2 ) = x^/2s following from (13. 8p . The signs of the complex 
conjugate values y/x~i, y/x^ must be chosen in such a way that the expressions (13.351) . (13.351) satisfy 
one of the equations (13.21) . (13.31) (then the other two hold automatically). 

Thus, all phase variables are algebraically expressed in terms of two auxiliary variables x, £; 
the domain of the latter depends on the constants of the first integrals. 
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4 Singularities of the induced symplectic structure 

The Hamiltonian structure of the system (11 .11) is provided by the Poisson brackets on R 9 (cj, a, f3). 
In notation (12 .41) these brackets are [2] 

{Mj, M k } = E jk iM h {Mj, a k } = e jM oi h {Mj, (3 k ) = e jH p h ^ . 

{a j; a fc } = {<x,, /3 fe } = {f3j,f3 k } = 0. 

Being restricted to P 6 they correspond to Kirillov's symplectic form A 6 A 2 (P 6 ). Recall the 
following well-known facts. Let N C P 6 be a submanifold defined in P 6 by two independent 
equations 

/i = 0, / 2 = (4.2) 

and let Xi , £ 2 be the Hamiltonian vector fields with the Hamilton functions f\ , f 2 . Then the 
span of £i,3l2 at eac h point £ G iV is skew orthogonal to the tangent space TqN. Therefore the 
restriction of A to iV is non-degenerate at the point ( if and only if 

{/i,/ 2 }(C) = A(X 1 ,X 2 )(C)^0. (4.3) 

Let us calculate the bracket {Pi, P 2 } of the functions fll . 13j) . As Pi is pure imaginary denote 

R = -{R 1 ,R 2 }. (4.4) 

1 

The change of variables (11.71) is linear with constant coefficients, so the rules (14. ip are easily 
transformed for new coordinates. Omitting technical details we present P in the form 

P = + F,f - P 4 ], (4.5) 

w 1 w 2 w 3 

where 

Pi = Z X Z 2 W 3 + X 2 ^iWi + X 1 Z 2 W 2 , 



P 2 = W1W2W3 + Z 2 Wi + Z 1 W 2 , 

P 3 = 2wfw2 + x 2 w\ + X\Vj\, 

2„„2|7„„2„„2 1 „ „„2 1 „ „„2> 



(4.6) 



P 4 = 4wfwl[(wfw% + X 2 w\ + XiW^Ws + X 2 Z 1 W 1 + XiZ 2 W 2 ]. 

The variables y±, y 2 have been eliminated as the solution of the system (11.121) 

WlU7 3 (x 2 Wi + Z 2 W 3 ) + Pi W 2 Wz(XiW 2 + ZiW 3 ) + Pi _s 

V\ — ) 2/2 = • (4.7) 

WiW 2 W 3 WiW 2 W 3 

Fixing the values h, k, g, s, r of the functions fll .3D . ( 12. lip . ( 12. 12ft we have 

w 3 F 2 + P 3 = 2w 1 w 2 h - 2{x x x 2 + z x z 2 - r), P 4 = Aw{w 2 2 w 3 (t - k). (4.8) 
Then from (T2TT41) . (IQl . (ET3I) . fCTBjl . (1331) we obtain 

^3P2 + P3 = f(^- S ), ^4 = -^, Pl=^3, ^2 = ^- ( 4 - 9 ) 



Finally, the expression (14.51) takes the form 



R= - 

s 



4 1 2 \ 2 1 P r 
s - {p - t)s z + 



(4.10) 



10 



Obviously, R is a first integral of the induced system, and the equation 

4 4 

s 4 - (p 2 - r)s 2 + = (4.11) 

defines the set of the integral constants s, r such that on the corresponding integral manifolds the 
2- form induced on D by the symplectic structure A degenerates. As the functions (12. lip . (12.121) 
have an algebraic structure, the set (14. lip has codimension 1 in O. Thus, \\q is almost everywhere 
non-degenerate. 

To locate the set ( 14.1ip on the surface ( 12. 141) in R 3 (/i, k,g) notice that, in virtue of ( 12. 141) . 

» (A, t, 9 ) x |(A, », 9) = p - - T y + £^1] <Iz£ ^, J,). 
Therefore, the equation (14. lip defines cuspidal edges of the surface ( I2.14p . 



5 Bifurcation diagram and the existence of motions 

Introduce the integral map J of the dynamical system on O 

J(C) = (S(C),T(C))eR 2 , (gD. (5.1) 

The bifurcation diagram E( J) of this map is, by definition, the set of pairs (s, r) over which J is 
not locally trivial. In our case all integral manifolds ( 13.11) are compact. Therefore, E( J) is the set 
of critical values of J. 

Define the admissible region as the image of the map (15.11) . i.e., the set of values (s, r) such 
that the integral manifold (13 .ip is not empty. Obviously, dJ(0) C S(J). 

Due to the above results the existence of a trajectory on O with given s, t is equivalent to 
the existence of a point on the (x, £)-plane for which the values (I3.26p - (l3.28p . (13.351) . (I3.35P satisfy 
(13. 7p . It easily follows from (I3.25p . (13.281) . (13. 8p that the conditions (13.71) . in turn, are equivalent 
to the system of inequalities 

$x(x,O^0, $ 2 (x,O^0, (5.2) 
^(x,O^0, # 2 (x,£)^0 (5.3) 

considered in the quadrant x ^ 0, s£ ^ 0. 

Theorem 3. The bifurcation diagram S(J) consists of the following subsets of the (s,r)-plane: 



1°) 


T = 


(a + b) 2 , s G [-a, 
(a — b) 2 , s G [—a, 


0) U [6,+oo); 


2°) 


T = 


-b] U (0,+oa); 


3°) 


S = 


—a, r ^ (a — b) 2 ; 




4°) 


S = 


-6, r ^ (a - b) 2 - 




5°) 


S = 


b, t<: (a + b) 2 ; 
a, t < (a + b) 2 ; 




6°) 


S = 




7°) 


T = 


0, s G (0,+oo); 




8°) 


T = 


(Va 2 - s 2 + Vb 2 - 


-s 2 ) 2 , se [-6,0); 


9°) 


T = 


(Va 2 - s 2 - Vb 2 - 


-^ 2 ) 2 , se(0,b}; 


io°; 


)r = 


-- -{Vs 2 -a 2 - Vi 
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Theorem 4. The solutions of the system (11.11) under the conditions (11.121) exist iff the constants 
of the first integrals (12. lip , f 1 2 . X 2 j) satisfy one of the following 

1°) - a ^ s < -b, t }t (a-b) 2 ; 

2°) -b < s < 0, r ^ ( y/a 2 - s 2 + ^fe 2 - s 2 ) 2 ; 

3°) < s ^ b, r <: (^a 2 -s 2 - v/6 2 - s 2 ) 2 ; 

4°) 6 < s < a, r ^ (a + b) 2 ; 

5°) s ^ a, -(^2-62 _ V / S 2 _ a 2)2 ^ r ^ ( a + fe )2 

The complete proof of these statements is purely technical (see [10]) and contains the scrupu- 
lous analysis of the regions on the (x, £)-plane defined by (15.21) . (15.31) and the cases of their 
structural transformations. Another approach is suggested in [TT] for the pair (S,H). It is based 
on the classification of the trajectories in D satisfying the condition rank (H x K x G) < 2. 









T 1 

(a+b) 2 








( 












■b) 2 




/ v 

\ 


s 




-a -b 






b a 


K 

r 



Figure 1: The admissible region in the (s, r)-plane. 



The admissible region is shaded in Fig.l. The dense lines and curves represent the equations 
of the bifurcation diagram, the dashed curve illustrates the equation (14.111) . i.e., the first integrals 
constants such that the symplectic structure is degenerate on the corresponding integral manifolds. 

6 Separation of variables 

In the sequel we suppose that to ^ 0. In fact, if r = 0, then from (I2.14p we obtain the relation 
(2g —p 2 h) 2 = r 4 k characteristic for the critical manifold 9T [9]. The equations of motion on 9T 
were explicitly integrated in [8]. If a — 0, then the equations (I2.14p yield one of the relations 
(11.111) . This case corresponds to the set of points (11.91) . At these points the manifold O fails to 
be smooth. The corresponding trajectories are the pendulum motions (11.101) . 

Considering the second equation (I3.2ip denote \i = \r 2 x\ — ryi\ = |r 2 x 2 — j - Then 

/i 2 = r £ 2 + ax 2 - to. (6.1) 

This equation defines the second-order surface M. in three-dimensional space R 3 (x, Each 
trajectory in O is in a natural way represented by a curve on Ji4. 
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Theorem 5. Supposing to 7^ 0, introduce the variables 



U 



V 



y/a(r - x 2 ) 
Then the equations of motion on O separate 

alU alV UdU 



VW) VW) 



0. 



y/a(r - x 2 ) ' 



VdV 



dt 



y/Q(U) yW) V^' 



Here 



Q( w ) = (w 2 - \)[aw 2 - As 2 x 2 )[(V^w + r) 2 - r 4 ]. 
Proof. Consider the local coordinates u, v on the surface Ai: 



uv + 1 



u — V 



uv 



U + V U + V U + V 

In addition to (13.191) note the following two identities for the constants introduced above, 



2\2 



a + 2r(p 2 ± r 2 ) = (r ± r # 
The polynomials (13.1 5p . (13.241) . (I3.36P become 

$1 = K^l(w) <Pl(^), $2 = K¥>2(«) ^2^), 
^1 = Klp^u) l[> 2 (v), #2 = Kll> 2 {u) i>l(v), 

Gi = K6 2 (u)d 1 (v), G 2 = k6 1 (u)6 2 (v), 

where k — l/(u + f ) 2 and 

ViH = v 7 ^ 1 + w 2 ) + 2(r + r 2 )w, <p 2 (w) = y/a(l + w 2 ) + 2(r - r 2 )w, 

^iH = 2s[(x + v^K - (X - v^)L ^O) = 2s[(x - v^)^ 2 - (x + v^)L 

0i(iu) = y/o(l - w 2 ) + Ast/tw, 6 2 (w) = y/a(l - w 2 ) - As^/tw. 

Then from f!3.26[) — (13. 28j) we find the expressions for the configuration variables 



Xl 



X-2 



2st 




+ <p 2 (u)tpi(v) 




_ v^i(«)*i(«) 


- V 2(u)9 2 (v) 


2st 


y/<Pi(u)<p 2 (v) 


- V V2(«)yi(«) 






+ y/e 2 (u)e 2 (v) 



-1 2 



2/1 = 2s 



2/2 = 2s 



-1 2 

y/<Pi(u)(p 2 (v) - y/ip2(u)ipi(v) -Aa(uv-lf 



1 



2l 

22 



2r(u + 1>) 
2r(ii + v) 



y/eMe^v) + y/e 2 (u)e 2 (v) 

V / ^lC^VlM+ \fV2{u)ip 2 {v 

y/(pi(u)ipi(v) - y/<p 2 (u)<p 2 (v 
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Hence, in particular, 



a/2sT y/<pi(u)<p 2 (v) + y/<P2(u)<Pi(v) 

r y /9 1 (u)9 1 (v) - y /6 2 {u)6 2 {v) 
y/2sr y/ipi(u)ip 2 (v) - y/<P2(u)<pi(v) 



(6.10) 



r y/9 1 {u)9 l {v) + ^/9 2 {u)9 2 {v) ' 
Here the arbitrary choice of sign is provided by the algebraic value \ / 2sr. From (I3.35j) we have 

45 U + V ' (6.11) 

>X\W 2 



X 2 Wi 



4s 



u + v 



Substitute (16. 10ft into ( 16. lip to obtain the expressions for the variables defining the equatorial 
components of the angular velocity 



ir 



Wi 



w 2 



VUu)9i(v) + v /9 l (u)9 2 (v) ^e l (u)9 1 (v) + y/0 2 (u)0 2 {v) 



As\/2sr 
ir 


(u + v) 




v) - ^W^iW 




VOiiuMv) - v /9 1 (u)9 2 (v) 


y/9i{u)e x (v) - ^e 2 (u)9 2 (v) 


W2sr (u + v) 


y/<pi(u)<p 2 (v) + \/<p 2 {u)(fi(v) 





The axial component is found from (I3.35p . 

i y/9 1 (u)8 2 (u)(pi(v)ip2(v) - ^ ip 1 {u)^> 2 {u)9 l (v)9 2 (v) 



w 3 



2V2 



srcr 



(u + v)(uv — 1) 



Thus we have expressed all phase variables in terms of u,v. 

To obtain the differential equations for u, v consider the following variables 



si 



x 2 + z 2 + r 2 
2x : 



s 2 



The derivatives, in virtue of the system ( 11. 8ft . are 

,.2 



x 2 + z 2 — r 2 
2x ' 



M - + ^X^l^ - a?2Wl), S 2 = |-3(2l - ^X^l^ + x 2 Wi). 



On the other hand, from (16.141) . (16. 5ft we have 

Ja(uv + 1) + (r + r 2 )(u + v) Ja(uv + 1) + (r - r 2 )(u + v) 
s i = 7T~F^ \ , -5 2 = 



2^[u - v) 



whence 



dsi 



dsi 



2y/¥{u - v) 



du 2y/r(u — f) 2 ' dv 2y/r(u — v ) 



ds 2 



(f 2 (v) 



ds? 



<p 2 (u) 



du 2y/r(u — v) 2 ' dv 2\[t{u — v) 1 



(6.12) 



(6.13) 



(6.14) 



(6.15) 
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Therefore, 

du 2Jt(u — v) 2 ds\ ,ds 2l 

dt ipi(u)ip 2 (v) - ip 2 (u)ipi(v) dt dt 

dv 2y/r(u — v) 2 dsi , ,ds 2 , 

~dt = Vl {u)^ 2 {v) - cp 2 (u)ipi(v) [(p2 ^~dt ~ ^^"rfp 

Substitute the values ( 16.9p ~ 116.lip into f ]6 . 1 5[) and the resulting expressions into ( 16. 16H . We obtain 



, .du y/<Pi(u)<p 2 (u) 9ju) e 2 (u) dv y/vi(v)<p 2 (v) 9jv) 6 2 (v) 

1{U ' V) H = • SM Tt = ' ( ' 

where 

. (u-v)(l — uv) ( 1\ ( 1 
f(u,v) = i ^ '- =(v + -)-(u + - 

uv V v / V u 

The variables ( 16. 2ft in terms of u, v have the form 

U = Uu + -), V = Uv + -). (6.18) 



2V u) 2V v 



To be definite choose u = U — \/U 2 — 1, v = V — \/V 2 — 1. Then the equations ( 16.171) yield 

(U-V)^ = ^J=yfQ(U), (U-V^^^^VW) (6-19) 
at \J2sto at y/Zsra 

with the polynomial ( 16.41) of degree 6. This system is obviously equivalent to the system ( 16.31) : 
the latter have the standard form for hyperelliptic quadratures. □ 

To reveal the connection of Theorem [5] with the bifurcation diagram suppose that the poly- 
nomial (16. 4p has a multiple root. The resultant of Q(w) and Q'(w) is 

2 3W(a 2 - & 2 ) WVxV - s 2 ) 2 (b 2 - s 2 ) 2 . (6.20) 

According to (12.141) . s^0on£). The equation o = gives r = (a ± b) 2 . The case x = provides 
the values r = (\/a 2 — s 2 ± a/& 2 — s 2 ) 2 . Thus, the bifurcation diagram of J = S x T belongs 
to the discriminant set of the polynomial (16. 4p . This fact is typical for systems with separating 
variables. 

Note that all roots of the polynomial ( 16. 4p are explicitly expressed in terms of the integral 
constants along with the roots of the polynomial 

<p 1 {w)ip 2 {w)6 1 {w)9 2 {w) (6.21) 

defining the solutions of the system ( 16. 17ft . The variables introduced in this section are real 
only in the case r > and o > 0. For other combinations of signs the definition (16. 5p of 
the local coordinates u, v needs obvious modifications. Then the intervals of oscillations for the 
variables U, V and u, v are easily determined for any given values of s, r or, more exactly, for each 
connected component of R 2 \S(J). The obtained differential equations and the explicit formulae 
for the phase variables in terms of u, v provide the possibility to investigate the phase topology 
of the case considered and to calculate analytically the numerical invariants of the corresponding 
Liouville foliation. 
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